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Phenyl-disubstituted polyacetylene (PDPA) is an organic semiconductor which has been studied 
during the last years for its efficient photo-luminescence. In contrast, the molecular geometry, 
providing the basis for the electronic and optical properties, has been hardly investigated. In 
this paper, we apply a density-functional-theory based molecular-dynamics approach to reveal the 
molecular structure of PDPA in detail. We find that oligomers of this material are limited in length, 
11^ . being stable only up to eight repeat units, while the polymer is energetically unfavorable. These 

facts, which are in excellent agreement with experimental findings, are explained through a detailed 
r- { I analysis of the bond lengths. A consequence of the latter is the appearance of pronounced torsion 

angles of the phenyl rings with respect to the plane of the polyene backbone, ranging from 55° up to 
95° . We point out that such large torsion angles do not destroy the conjugation of the n electrons 
from the backbone to the side phenyl rings, as is evident from the electronic charge density. 



PACS numbers: 71.15.Mb, 71.20.Rv, 31.15.Qg 



INTRODUCTION 



. In the emerging field of nanotechnology conjugated polymers are playing an important role as materials of choice 

jrt for organic light emitting diodes (OLEDs), (author?) ^'^ organic field effect transistors (OFETs), (author?) - organic 

Cj lasers, (author?) '''^ photocells, (author?) '^ etc. Their application in such devices requires the search and the inves- 

i' \ tigation of new materials with interesting optical and electronic properties. PA which is the simplest and the most 

'^ ■ widely studied conjugated polymer, exhibits very weak photo-luminescence(author?)— ili^i (PL), thus ruling it out as 

^ I a candidate for opto-electronic applications. However, in the last decade, phenyl-disubstituted polyacetylene (PDPA) 

w : derivatives, obtained by replacing hydrogen atoms of irans-polyacetylene (PA) by phenyl rings or their derivatives, 

were demonstrated as luminescent materials with high quantum efficiency, which also exhibited stimulated emission 

in thin film form, (author ?)—iiS Therefore, one anticipates that they will be useful in creating light emitting diodes 

and polymeric lasers. (author?) ^1 

^ , To understand the physics behind the fluorescence of these materials, a series of experimental and theoretical 

investigations have been performed, (author?)— J^J^J^J^J^J^Jiii^J^i^i^i^ From calculations based on a Pariser- 

Parr-Pople (PPP) model Hamiltonian it was argued that this phenomenon is due to reverse excited-state ordering 

(compared to PA) of the one-photon allowed IBu and two-photon allowed 2Ag states, caused by reduced electron 

correlation effects due to the presence of phenyl rings. (author?)— J^Ji In PA, the 2Ag state occurs below the IBu 

level, as a consequence of which the optically pumped 1_B„ state decays rapidly to the 2Ag state, and according to 

^^ ' dipole transition rules, radiative transition between 2Ag and the lAg states is forbidden. Hence it was concluded that 

in contrast to PA, PDPA shows strong PL as E{2Ag) > E{lBu), where li?„ is strongly dipole coupled to the ground 

state. 

While most of the investigations concentrated on the optical properties of PDPAs their structural properties have 
hardly been investigated. A complete understanding of the electronic and hence optical properties, is, however, 
5— i impossible without knowing the respective ground state geometries. This leads to another puzzling issue, namely 

the reason behind the experimentally observed short conjugation length of these materials, which typically consist of 
seven repeat units in PDPA thin films, (author?) ^^ and only five in solution. (author?) ^^ Based upon an SSH model 
calculations, it was predicted in an earlier work that for an infinite chain, introduction of phenyl rings will lead to a 
reduction of bond alternation along the backbone. (author?) ^^ 

The aim of this work is to throw light on the formation of PDPAs by investigating their energetics as a function of 
oligomer length. To this extent we carry out first-principles calculations based on density functional theory (DFT). We 
thereby focus on the geometry, including a detailed study of bond lengths and alternations, as well as the orientation 
of the phenyl rings with respect to the polyene backbone. We will show that polymerization of PDPSs is hampered 
by single-bond breaking, and oligomers are only stable up to eight repeat units. Moreover, an analysis of band gaps 
and charge densities allows for an insight into the lowest optical transitions. 




Figure 1: (color online) Repeat units of poly acetylene and phenyl-disubstituted polyacetylene. Large and small spheres denote 
carbon and hydrogen atoms, respectively. Ts represent single bonds and rd double bond lengths. 

II. PDPA AND ITS OLIGOMERS 

The crystal structure of PA is of herringbone type, where the polymer chains are coupled by the weak van der 
Waals interactions. Although the 3D environment can have a large impact on the optical properties on such sys- 
tems, (author?) ^"^'^^ it can be considered as quasi one-dimensional for the current investigation, where we study 
the cohesive properties as a function of chain length. Hence we model PA, and consequently PDPA, as strictly 
one-dimensional zigzag chains. 

The PDPAs with the formula unit C14H10 are obtained by replacing the hydrogen atoms of PA (C2H2) by phenyl 
groups, leading to the structure presented in Fig. [TJ The backbone polyene chain is running along the x axis 
(longitudinal direction) while the phenyl rings are oriented along the y axis (transverse direction). The single and 
double bonds are labeled as r^ and r^, respectively, while the bond alternation is defined as Ar = rg — rd- a and j3 
denote the angles between the carbon atoms of the backbone and the phenyl ring, 7 is the angle between the carbon 
atoms of the backbone, and 5 is the torsion angle of the phenyl ring with respect to the backbone plane. The PDPA 
oligomers are terminated with one hydrogen atom at either end, saturating the dangling bonds. The formula units of 
PDPAs and PAs are Ci4„Hion+2 and C2nH2„+2, respectively, with n being a positive integer representing the number 
of repeat units. We adopt the notation PDPA-ra. and correspondingly PA-n for the oligomers of PA. 

III. COMPUTATIONAL DETAILS 

In the present work, infinite PDPA as well as its oligomers with two to six building blocks are investigated by 
means of density functional theory employing the supercell approach. A vacuum region of around 8 A is adopted 
along the y and z directions for the polymers, and also in the third dimension for the oligomers. This size is 
considered sufficient to avoid interaction between the translational images. All the geometries are optimized, where 
we assume the structures to be relaxed when the forces on individual atoms are below 1 meV/A. To this extent, 
the projector augmented wave (PAW) method is utilized as implemented in the original PAW code. (author?) — 
It makes use of Car-Parinello molecular dynamics allowing for geometry optimization in an accurate and efficient 
manner. Polynomial-type pseudopotentials as implemented in the code are used, with a plane wave cutoff of 30 Ry. 
A Brillouin zone sampling of 16 x 1 x 1 k points is taken for the calculations of polymers, while the finite systems 
naturally requires one k point only. Exchange-correlation effects are treated by the generalized gradient approximation 
(GGA) in the Perdew-Burke-Ernzerhof (PBE) flavor, (author?)— We further discuss the issue of our choice of the 
exchange-correlation functional at the end of section IIVAI 

The polymerization energy Epoiy is calculated using the formula 

^poly — ^polymer ^monomer V-tJ 



where Epoiymer is the total energy per unit cell of the infinite polymer and Emonomer is that of diphenyl acetylene, 
C14H10, in its optimized geometry. The choice of C14H10 as the monomer is based upon the fact that it has the same 
chemical formula as of the repeat unit shown in Fig. [1] Moreover, in experimental situation also PDPA derivatives 
are usually synthesized using diphenyl acetylene monomer and the TaC15-Bu4Sn catalyst. (author?)-- Therefore, we 
believe that our choice of using diphenyl acetylene as a monomer is on a strong footing from a chemical point of view. 



Similarly, we define an oligomerization energy, EoUgo, as 

1 



n 
where EoUgomer is the total energy of the respective oligomer and Eh is that of the hydrogen atom. 



(2) 



IV. RESULTS 



A. Structural Properties 



To evaluate the stability of PDPA we calculate the polymerization energy as a function of the lattice parameter a, 
which is displayed in Fig. [5] (bottom panel). The most striking result is that the values are positive over the whole 
range, indicating that the system is not bound at all. The energy monotonically decreases with increasing lattice 
parameter up to a « 3.1 A followed by a double kink which is explained by a change in the bonding characteristics. On 
analyzing the single and double bond lengths (top panel) we see that first both of them increase due to steric hindrance 
between the phenyl rings, until the double bond reaches a value of « 1.45 A. With even larger lattice spacing, the 
single bond length keeps increasing while the double bond length goes down, because any further increase would 
correspond to a single bond, which is impossible given the chemical structure. Once a reaches 3.115 A, r^ becomes 
smaller than 1.34 A, while r^ gets larger than 1.54 A, values typical for C=C and C-C bonds, respectively. This finally 
leads to the breaking of the single bond and the subsequent transformation of the double bond into a triple bond, 
converting the system into individual diphenyl acetylene molecules as indicated by the fact that the binding energy 
is approaching zero. 
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Figure 2: Variation of the single bond length, r^, double bond length, 
PDPA (bottom) as a function of the lattice parameter a. 



r-d 



(top) and polymerization energy, Epoiy, of infinite 



For comparison, we perform the same calculation for PA. Kahlert and co-workers (author?)— measured the lattice 
parameter as well as C=C and C-C bond lengths to be 2.455 A, 1.36 A, and 1.45 A, respectively, while Yannono 
et al. (author?) ■^^ reported on the basis of NMR results double and single bond lengths of 1.36 A and 1.44 A, 
respectively. Our results, shown in Fig. [3l demonstrate, as expected, that the system is bound, with the theoretical 
lattice parameter being 2.475 A, and the corresponding values of r^ = 1.42 A, r^ = 1.38 A. Both, theory and 
experiment exhibit nonzero bond alternation, which is fully consistent with Peierl's theorem, (author?) '^^ stating that 
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Figure 3: Variation of the single bond length, r^, double bond length, r^, (top) and polymerization energy, Epoiy, of polyacety- 
lene, as a function of the lattice parameter a. 



one-dimensional metals undergo a lattice distortion to minimize their energies, thus becoming band-gap materials. 
However, the magnitude of bond alternation, Ar = 0.04 A, found by our calculations is roughly a factor of two 
smaller than the experimental value, (author?)— *^ This underestimation of bond alternation by DFT calculations 
is a well-known problem, discussed earlier by various authors. (author?) "^^ 

Realizing that the PDPA polymer is not stable, the question arises what limits the length of its oligomers. To this 
extent we study the energetics as a function of molecular size, considering two to six repeat units. For each of these 
structures, the geometry is optimized and the results are summarized in Table ID We note that all the oligomers are 
found to be invariant under inversion symmetry, exhibiting point group Ci. 

Their oligomerization energies, EoUgo, as defined earlier, are displayed for these representatives in Fig. |4l ex- 
hibiting a decrease with increasing molecular size. To highlight the instability with larger oligomer length, we have 
plotted the energy as a function of inverse length, i.e., 1/n in the inset of the figure. Here, extrapolation clearly 
indicates that the formation energy will become positive for n > 8. This is in excellent agreement with experimental 
findings, (author?)— where the longest oligomers have been estimated to consist of not more than seven repeat units. 

It will be interesting to perform these calculations also for PDPA-7, however, given the computational constraints, 
at present it is not feasible. 

Analyzing the structural parameters, we observe a clear trend for the bond distances. Single as well as double bond 
lengths not only increase with the oligomer length, but also for a given oligomer they increase going from the outer 
PDPA repeat units toward the innermost ones. For example, PDPA-6 rj=1.39 A and r^=1.44 A differ by as much as 
0.05 A while Ts increases from 1.52 A to 1.56 A at the same time. The latter value signalizes the bond breaking for 
larger oligomers as found for the polymer before. In contrast, for the PA oligomers the corresponding values all lie in 
the range of the polymer bond lengths exhibiting only small variations. The different bond lengths are reflected in the 
behavior of the bond alternation along the backbone. For example, Ar for the innermost bonds in PDPA-6 is 0.12 A, 
while it is only half as large in PA-6. These results are in clear contrast to previous model calculations. (author?) ^~ 

The instability of longer oligomers is also reflected by the orientation of the phenyl rings with respect to the 
backbone, what is evident from Fig. [S] showing the relaxed geometry of PDPA-6. While at the innermost repeat units 
the phenyl groups take a nearly standing position with respect to the x axis, the outer ones are leaning toward the 
outside. This is demonstrated by the angle a {(3), dramatically decreasing (increasing) from the left to the right. 

Another interesting aspect concerns the torsion angle, 6, of the phenyl ring with respect to the polyene plane. In 
contrast to previous work assuming it close to 30°, (author?) ^^ we predict much larger values, ranging from 56° to 
more than 90°. This appears counter- intuitive at the first sight, because one would expect that such large torsion 
angles would destroy the conjugation between the backbone and the phenyl rings. This fact is, however, required 
by the small spacing between the rings provided by the polyene backbone. Only large torsion angles allow for an 
accommodation of the phenyl groups substituting the hydrogen atoms. A similar behavior is seen in rubrene, in which 
equally large torsion angles for the side phenyl rings have been reported. (author?) ^2,33 

Before closing the discussion of the structural properties of PDPA's, we would like to briefly comment on the 



Table I: Structural parameters for the oligomers of PDPA and PA, respectively. For definitions of the angles see Figure [T] Bond 
lengths Td and r^ as well as bond alternations Ar are given in A. The superscripts denote the index of the repeat unit. 





PDPA-2 


PDPA-3 


PDPA-4 


PDPA-5 


PDPA-6 


PA-2 


PA-3 


PA-4 


PA-5 


PA-6 


ri 


1.37 


1.38 


1.38 


1.38 


1.39 


1.35 


1.35 


1.35 


1.35 


1.35 


2 
Td 


1.37 


1.40 


1.41 


1.42 


1.43 


1.35 


1.36 


1.37 


1.37 


1.37 


ri 




1.38 


1.41 


1.43 


1.44 




1.35 


1.37 


1.37 


1.37 


ri 






1.38 


1.42 


1.44 






1.35 


1.37 


1.37 


rl 








1.38 


1.43 








1.35 


1.37 


rl 










1.39 










1.35 


rl 


1.48 


1.49 


1.50 


1.51 


1.52 


1.45 


1.44 


1.44 


1.44 


1.44 


2 




1.49 


1.51 


1.53 


1.55 




1.44 


1.44 


1.43 


1.43 


rl 






1.50 


1.53 


1.56 






1.44 


1.43 


1.43 


rt 








1.51 


1.55 








1.44 


1.43 


rl 










1.52 










1.44 


Ar' 


0.11 


0.11 


0.12 


0.13 


0.13 


0.10 


0.09 


0.09 


0.09 


0.09 


Ar^ 




0.09 


0.10 


0.11 


0.12 




0.08 


0.07 


0.06 


0.06 


Ar^ 






0.09 


0.10 


0.12 






0.07 


0.06 


0.06 


Ar* 








0.09 


0.11 








0.07 


0.06 


Ar^ 










0.09 










0.07 


a' 


134.0° 


138.8° 


141.6° 


143.3° 


143.4° 


121.2° 


121.1° 


121.2° 


121.1° 


121.3° 


2 

a 


113.2° 


120.8° 


124.7° 


126.7° 


127.6° 


119.5° 


118.7° 


118.8° 


118.8° 


118.7° 


a^ 




108.3° 


115.4° 


119.0° 


120.1° 




119.1° 


118.5° 


118.5° 


118.4° 


a" 






105.6° 


112.6° 


115.3° 






119.0° 


118.4° 


118.4° 


a' 








104.0° 


111.0° 








119.0 ° 


118.3° 


a« 










103.6° 










119.0° 


P' 


112.2° 


109.2° 


107.2° 


105.9° 


105.9° 


117.1° 


117.2° 


117.1° 


117.1° 


117.1° 


P' 


119.8° 


112.7° 


108.9° 


106.5° 


105.7° 


116.4° 


116.8° 


117.0° 


116.9° 


117.1° 


P' 




125.3° 


117.2° 


112.9° 


110.8° 




116.3° 


117.0° 


117.2° 


117.3° 


P' 






128.20° 


120.1° 


115.8° 






116.4° 


117.0° 


117.1° 


P" 








129.8° 


121.6° 








116.5° 


117.2° 


P' 










130.3° 










116.4° 


l' 


113.8° 


112.0° 


111.2° 


110.8° 


110.7° 


121.7° 


121.7° 


121.7° 


121.7° 


121.6° 


7^ 


127.0° 


126.5° 


126.4° 


126.8° 


126.7° 


124.1° 


124.5° 


124.2° 


124.3° 


124.2° 


7^ 




126.4° 


127.4° 


128.1° 


129.1° 




124.6° 


124.5° 


124.3° 


124.3° 


7^ 






136.4° 


127.3° 


128.9° 






124.6° 


124.6° 


124.5° 


7^ 








126.2° 


127.4° 








124.5° 


124.5° 


7*^ 










126.1° 










124.6° 


&' 


56.3° 


65.7° 


72.4° 


80.6° 


78.0° 












5^ 


63.4° 


77.1° 


85.9° 


92.1° 


86.9° 












5^ 




68.1° 


78.0° 


84.1° 


81.9° 












5^ 






71.3° 


77.8° 


79.2° 












S^ 








74.0° 


75.8° 












5" 










71.2° 













choice of the exchange-correlation functional used in these calculations, vis-a-vis some of the hybrid functionals which 
have become popular these days, (author?)— These hybrid functionals contain a certain percentage (20%-40%) of 
exact HF exchange and give quite good results for conjugated polymers in several cases, (author?) ■^'' but there is 
no known value of the exchange percentage which will give correct results for all the properties. To investiage as to 
whehter by using a hybrid functional one would obtain substantially different results than our GGA-PBE approach, 
we performed these structural calculations on PDPA-2 using the PBEO hybrid functional(author?) ^^ as implemented 
in the ESPRESSO program(author?) ^^ with 25% HF exchange. The calculation time for the PDPA-2 molecule with 
the PBEO functional increases 400 times, as compared to the GGA-PBE method, thereby making such calculations 
impossible for longer oligomers. The structural parameters and the oligomerization energy of PDPA-2 obtained using 
these two functionals are presented in Table |lll and from the results it is obvious that most of them agree almost 
exactly with each other. The only notable differences are: (a) 3.5% disagreement in the value of 5' , and (b) about 9% 
difference in EoUgo, which, for the present calculations will not alter our conclusions in any significant manner. We 
believe that similar trends will also hold for longer oligomers, and, therefore, the choice of the GGA-PBE approach 
in our case is fairly adequate. 
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Figure 4: Oligomerization energy of PDPA as a function of oligomer length. The inset presents the same data with respect to 
the inverse number of repeat units. The full line represents a cubic-spline fit to the calculated values, indicating that structures 
longer than eight repeat units will be unstable. 




Figure 5: (color online) Relaxed geometry of PDPA-6. The angles are exemplary labeled for the third repeat unit. Note the 
inversion symmetry of the system. 



B. Optical Absorption 



Having discussed the structural properties, we now have a look at their impact on the electronic structure. To 
this extent, we compare the gap between the highest occupied molecular orbital (HOMO) and the lowest unoccupied 
molecular orbital (LUMO) of oligo-PDPA's and polyenes, presented in Table iHll As expected, in both the systems, 
the energy gap decreases with the increase in molecular length. However, the gap is always found to be smaller 
in case of PDPA-n compared to PA-n, which is in agreement with experimental results(author?) ^'^° and earlier 
calculations. (author?) ^^ Note that the presented values are Kohn-Sham gaps with the PBE exchange-correlation 
functional, (author?) ^^ which are well known to be underestimated by roughly a factor of two. 

In order to ensure that the Kohn-Sham gaps actually correspond to the optical gaps, we calculated both the 

(2) (2) 

longitudinal and the tranverse components of the imaginary part of dielectric constant tensor, exx{oj) and Cyy (w), 
respectively {x being the conjugation direction) for PDPA-2, and the results are presented in Fig. [HI From the results it 



A2), 



J2) 



is obvious that: (a) the first peaks in both ei,^; (w) and eyy{uj) occur at 2.88 eV, precisely the Kohn-Sham gap of PDPA- 
2 reported in Table iHll and (b) significant intensity of the first peak of e^j, (w), indicating a non-trivial transverse 
polarization component in the linear absorption at the optical gap, consistent with the earlier works, (author?)— 



Table II: Structural parameters and oligomerization energy of PDPA-2 calculated by GGA-PBE and the PBEO (hybrid func- 
tional) approaches. For definitions of the angles see Figure [T] All the lengths are in A, and the energies in niilliHartrees. 
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Table III: Kohn-Sham band gaps for various oligomers of PDPA and PA. 
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In order to account for the excitonic effects in the optical absorption, more sophisticated calculations including the 
effects of electron-hole interaction are needed, which are outside the scope of the present work. Nevertheless, we 
expect that in parallel to the gaps, the exciton binding energies will also decrease with the increasing oligomer 
length, (author?)^^^ 

Next we examine the the presence of the significant transverse polarization in the HOMO-LUMO absorp- 
tion(author?) ^^'^^ discussed above, from the point of view of the charge densities of these frontier orbitals. In 
light of the large torsion angles of the phenyl rings this observation needs to be better understood. For this purpose 
we present the wave functions of the HOMO and LUMO in Fig. [T] While the HOMO has significant charge density, 
both on the backbone as well as on the phenyl rings, the LUMO charge is much more confined to the backbone. This 
can be better understood when taking the difference between the two corresponding densities. Thus, HOMO— >LUMO 
transitions in the optical spectra should involve significant charge transfer from the side groups to the backbone atoms, 
giving rise to transverse polarization. Moreover, the fact that the HOMO has substantial charge density on the phenyl 
rings also indicates that the lar '" ^- -■- 1-- j- -~^ j--^ ^i- ~ -i--^ — — • ^■--- ^jetween the backbone and 




.(2) 



Figure 6: (color online) Longitudinal (e^J (oj)) and transverse components (e^y (ui)) of the imaginary part of the dielectric 



(2)/ 



constants of PDPA-2 plotted as a function of the energy of incident radiation. 
PBE(author?)— exchange-correlation functional. 



The calculations were performed using the 




Figure 7: (color online) HOMO and LUMO wave functions (isovalue of ±0.01), and corresponding charge density difference 
(isovalue of ±0.001 e) plot. Dark (light) regions denote positive (negative) charge. 



the phenyl rings. A similar trend was also seen by Petrenko et at for the case of rubrene. (author?)' 
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V. CONCLUSIONS 

In summary, we have demonstrated by total-energy calculations based on density functional theory that PDPA 
oligomers are only stable up to eight repeat units, while the corresponding polymer cannot be formed. This finding 
is related to the large size of the phenyl rings which need to be accommodated within the length of a polyene dimer 
(repeat unit of the backbone) confirming and explaining experimental observations. The torsion angles of the phenyl 
rings range from 55° to 95°, ruling out an earlier conjecture about their values being 30° only. At the same time we 
have shown that such large torsion angles do not destroy the conjugation of electrons from the backbone with those of 
the phenyl rings. This fact is supported by the electron distribution of the HOMO which exhibits substantial charge 
density on the phenyl rings. 

The energetics and the resulting structural properties of the the oligo-PDPAs have been compared in detail with the 
corresponding situation in the parent-polymer, polyacetylene. It has been revealed that due to phenyl substitution, 
the bond alternation increases while the band gap decreases. 
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